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Abstract 

Prehypertension is a risk factor for atlierosclerosis. We investigated alterations in plasma metabolites that are associated 
with prehypertension. A group of 53 individuals was identified who remained within the range of prehypertension during 
repeated measurements in a 3-year period. This group was compared with the control group of 53 normotensive subjects 
who were matched for age and gender. Metabolomic profiles were analyzed with UPLC-LTQ-Orbitrap mass spectrometry. 
The prehypertensive group showed higher levels of lysophosphatidylcholines (lysoPCs) containing C14:0, C16:1, C16:0, 
C18:2, C18:1, C18:0, C20:5, C20:4, C20:3, and C22:6, higher circulating Lp-PLA2 activity, oxidized LDL (ox-LDL), interleukin 6 
(IL-6), urinary 8-epi-PGF2a, and higher brachial-ankle pulse wave velocity (ba-PWV), before and after adjusting for BMI, WHR, 
smoking, alcohol consumption, serum lipid profiles, glucose, and insulin. LysoPC (16:0) was the most important plasma 
metabolite for evaluating the difference between control and prehypertensive groups, with a variable important in the 
projection (VIP) value of 17.173, and it showed a positive and independent association with DBP and SBP. In the 
prehypertensive group, the levels of lysoPC (16:0) positively and significantly correlated with ox-LDL, Lp-PLA2 activity, 8-epi- 
PGF2a, ba-PWV, and IL-6 before and after adjusting for confounding variables. Prehypertension-associated elevations in 
lysoPCs, Lp-PLA2 activity, ox-LDL, urinary 8-epi-PGF2ri, IL-6, and ba-PWV could indicate increased oxidative stress from Lp- 
PLA2-catalyzed PC hydrolysis during increased LDL oxidation, thereby enhancing proinflammation and arterial stiffness. 
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Introduction 

Hypertension is a risk factor for atherosclerosis and cardiovas- 
cular disease (CVD) [1-4], although the mechanisms by which 
hypertension is related to atherosclerosis are not clearly estab- 
lished. Several metabolomic studies have been published that 
investigate the metabolic effects of hypertension [5-7]. These 
studies reported abnormalities in gender-linked steroid patterns [5] 
or lipid metabolism [6,7]. Hypertension and its underlying 
pathophysiology may be present for years before clinical diagnosis, 
at which time irreversible pathology has already occurred. In 
2003, the Seventh Joint National Committee on Prevention, 
Detection, Evaluation, and Treatment of High Blood Pressure 
created a prehypertension category for persons with blood 
pressures ranging from 120-139 mm Hg (systolic) or from 80- 
89 mm Hg (diastolic). The committee strongly advocated lifestyle 
and behavioral modifications for individuals with prehypertension 
[8] . Prehypertension can precede hypertension and atherosclerosis 
for decades, and it is a condition that represents early CVD. 
Therefore, it is necessary to determine the role of prehypertension- 
associated alterations in circulating metabolic profiles. We 
performed metabolic profiling in a group of 53 individuals who 



remained within the range of prehypertension during repeated 
measurements in a 3-year period, and compared these with die 
metabolic profiles of age- and sex-matched normotensive controls 
in the same cohort. We also determined lipoprotein-associated 
phospholipase A2 (LP-PLA2) activity, oxidized LDL (ox-LDL), 
lipid peroxides, and brachial-ankle pulse wave velocities (ba- 
PWV). 

Materials and Methods 

Study subjects 

The 3-year prospective cohort study included 600 healthy 
subjects (30-65 years old) who underwent triennial medical 
evaluation from January 2008 to December 2011 through the 
National Health Insurance Corporation llsan Hospital in Goyang, 
Korea. Prehypertension was defined as blood pressure of 120- 
1 39 mm Hg (systolic) and 80-89 mm Hg (diastolic). The screening 
identified 53 nonobese individuals with normal fasting glucose (29 
men and 24 women) who remained within the range for 
prehypertension during repeated measurements in the 3-year 
period. From the same population, 53 normotensive control 
individuals matched for age and sex were recruited who blood 
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pressure has remained within the normal range (systoUc blood 
pressure <120 mm Hg and diastolic blood pressure <80 mm Hg) 
during repeated measurements in the 3-year period. The purpose 
of the study was carefully explained to all participants, and written 
consent was obtained prior to their participation. The Institutional 
Review Board of the NHIC-sponsored Ilsan Hospital and Yonsci 
University provided ethical approval of the study protocol, which 
was performed according to the Helsinki Declaration. 

Anthropometric parameters, blood pressure, blood 
collection, and dietary intake assessments 

Body weight and height of unclothed subjects without shoes 
were measured in the morning to calculate body mass index (BMI, 
kg/m^). Waist circumference was measured at the umbilical level 
with the subjects standing after normal expiration. Blood pressure 
(BP) was measured in the left arm of seated patients with an 
automatic BP monitor (TM-2654, A&D, Tokyo, Japan) after a 20- 
minute rest. After a 12-hour fasting period, venous blood 
specimens were collected in EDTA-treated and plain tubes, 
centrifuged to produce plasma or serum, and stored at — 70°C 
until analysis. The usual dietary intake of the study subjects was 
assessed using a semi-quantitative food-frequency rjuestionnaire 
and a 24-hour recall method [9]. Nutrient intake was dc-tc-rmined 
and calculated based on three-day food records using the 
Computer- Aided Nutritional Analysis Program (CAN-pro 2.0; 
Korean Nutrition Society, Seoul, Korea). 

Serum lipid profile and free fatty acids 

Fasting levels of total cholesterol and triglyceride were measured 
using commercially available kits and a Hitachi 7150 Autoanalyzer 

(Hitachi Ltd., Tokyo, Japan). ApoB-containing lipoproteins were 
precipitated with dextran-sulfate magnesium, and HDL-cholester- 
ol concentrations in the supernatants were measured enzymati- 
caUy. For subjects with serum triglyceride levels <400 mg/dL, 
LDL-cholesterol concentrations were estimated indirecdy using 
the Friedwald formula: LDL-cholesterol — Total-cholesterol - 
[HDL-cholesterol + (Triglycerides/5)] . For subjects with serum 
triglyceride levels a40() mg/dL, LDL-cholesterol concentrations 
were measured indirecdy. Free fatty acids were analyzed with a 
Hitachi 7150 Autoanalyzer (Hitachi Ltd, Tokyo, Japan). 

Fasting glucose, insulin, and homeostasis-model 
assessment of insulin resistance 

Fasting glucose levels were measured by the glucose-oxidase 
method with a Beckman Glucose Analyzer (Beckman Instruments, 
Ir\dne, CA, USA). Insulin levels were measured by radioimmu- 
noassay using commercial kits from Immuno Nuc:leo Corporation 
(StUlwater, MN, USA). Insulin resistance (IR) was calculated by 
the homeostasis-model assessment (HOMA) using the following 
equation: HOMA — IR = [Fasting insulin (nIU/mL) x Fasting 
glucose (mmol/L)]/22.5. 

Assessment of serum high-sensitivity C-reactive protein, 
lipoprotein-associated PLA2 activity, plasma 
malondialdehyde, and LDL particle size 

The concentrations of serum high-sensiti\"ity C-reactive protein 
(hs-CRP) were measured with an Express Plus TM auto-analyzer 
(Chiron Diagnostics Co., Walpole, NIA, USA) using a commer- 
cially available, high-sensitivity CRP-Latex(II) X2 kit (Seiken 
Laboratories Ltd., Tokyo, Japan). The activity of lipoprotein- 
associated phospholipase A2 (Lp-PLA2) was measured using a 
modification of a previously described, high-throughput radio- 
metric activity assay [10]. Plasma malondialdehyde (MDA) was 



measured from thiobarbituric acid-reactive substances using the 
TEARS Assay Kit (Zepto-Metrix Co., Buffalo, NY, USA). LDL 
particles were isolated by sequential flotation ultracentrifugation; 

the particle size distribution (1.019-1.063 g/mL) was examined by 
a pore-gradient lipoprotein system (CBS Scientific, CA, USA) on 
commercially available, non-denaturing, polyacrylamide slab gels 
containing a linear gradient of 2-16% acrylamide (Alamo Gels 
Inc., San Antonio, TX, USA). Latex bead (34 nm) conjugated 
standards of thyroglobulin (17 nm), apoferritin (12.2 nm), and 
catalase (10.4 nm) were used to estimate the relative migration 
rates of each band. Gels were scanned using a GS-800 Calibrated 
Imaging Densitometer (Bio-Rad, Graz, Austria). 

Plasma oxidized LDL, urinary 8-ep/-prostaglandi 
serum interleukin-6, interleukin-1 p, tumor necrosis 
factor-ot, and measurement of brachial-ankle pulse wave 
velocity 

Plasma ox-LDL was measured using an enzyme immunoassay 
(Mercodia, Uppsala, Sweden). The resulting color reaction was 
determined at 450 nm on a Wallac Victor^ multilabel counter 
(Perkin-Ekner Life Sciences, Turku, Finland). Urine was collected 
in polyethylene bottles containing 1 % butylated hydroxytoluene 
after a 12-hour fast. The botties were immediately covered with 
aluminum foil and stored at — 70°C until further analysis. The 
compound S-^jfJZ-PGFjn was measured using an enzyme immuno- 
assay (Bioxytech urinary 8-epi-VGF2oi Assay kit, OXIS Interna- 
tional Inc., Portiand, OR). Urinary creatinine levels were 
determined using the alkaline-picrate Jaffe) reaction. Serum 
interleukin (IL)- 1 fi. IL-6, and tumor nc'crosis factor (TNF)— a 
concentrations were measured using the Bio-Plex Reagent Kit and 
a Bio-Plex (Bio-Rad Laboratories, Hercules, CA, USA) according 
to the manufacturer's instructions. Brachial-ankle pulse wave 
velocity was measured using an automatic waveform analyzer 
(Model VP- 1000; Nippon Colin Ltd., Komaki, Japan) according to 
a previously published method [11]. The average ba-PWV from 
both left and right sides was used for analysis. 

Global (nontargeted) metabolic profiling of plasma 

Sample preparation and analysis. Prior to analysis, 
800 |xL of 80% acetonitrile was added to 100 |xL of plasma, 
mixed by vortexing, and centrifuged at 1 0,000 rpm for 5 minutes 
at 4°C. The supernatant was dried with Nj, dissolved in 10% 
methanol, mixed by vortexing, and centrifuged at 10,000 rpm for 
5 minutes at 4°C. The supernatant was transferred into a vial. 

Ultra performance liquid chromatography. The plasma 
extract samples (7 )j.L) were injected into an Acquity UPLC-BEH- 
C18 column (2.1 x 50 mm, 1.7 nm; Waters, Milford, MA) that was 
coupled in-line with a UPLC-LTQ_-Orbitrap XL (Thermo Fisher 
Scientific, USA). The injected samples were equilibrated with 
water containing 0.1% formic acid. Samples were eluted with an 
acetonitrile gradient containing 0. 1 % formic acid at a flow rate of 
0.35 mL/min for 20 minutes. Metabolites were separated by 
UPLC (Waters, Milford, MA), analyzed, and assigned by LTQ; 
Orbitrap-XL (Thermo Fisher Scientific, USA). The mass 
spectrometer was operated in ESI-positive mode. The spray 
voltage was 5 kV. The flow-rate nitrogen sheath gas and the 
auxiliary gas were 50 and 5 (arbitrary units). The capillary voltage 
(V), tube-lens voltage (V), and capillary temperature (°C) were kept 
constant at 35 V, 80 V, and 370°C. The Orbitrap data were 
collected in the range of m/z 50-1,000. For quality control, a 
mixture of four standard compounds (acetaminophen, sulfadi- 
methoxine, terfenadine, and reserpine) was injected every ten 
samples. The MS/MS spectra of metabolites were obtained by a 
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collision-energy ramp from 55-65 eV, and conducted with 
Xcalibur 2. 1 and MS Frontier software (Thermo Fisher Scientific, 
USA). 

Data processing and identification of metabolites. All 

MS data including retention times, ml z, and ion intensities were 
extracted hy SIEVE s()ftwar(; (Tliermo Fisher Scientific, USA) 
incorporated into the instrument, and the resulting MS data were 
assembled into a matrix. SIEVE parameters were set as follows: 
m/ z range 50-1,000; m/ z width 0.02; retention-time width 2.5; 
and m/ z tolerance 0.005. Metabolites were searched using the 
following databases: ChemSpider (www.chemspider.com), Human 
Metabolome (www.hmdb.ca), Lipid MAPS (v\'ww.lipidmaps.org), 
KEGG (www.genome.jp/kegg), and Ma.ssBank (www.massbank. 
jp). Selected metabolites were confirmed on the basis of retention 
times and mass spectra of standard samples. 

Statistical analyses 

Statistical analyses were performed using SPSS v. 21.0 (IBM 
SPSS Statistics 21, Chicago, IL, USA). Skewed variables were 
logarithmically-transformed for statistical analyses. For descriptive 
purposes, mean values are presented using untransformed values. 
Results are expressed as means ± standard error (SE). A two-tailed 
P-value of <0.05 was considered statistically significant. Differ- 
ences in clinical variables between the normotensive and 
prehypertensive groups at the three-year foUow-up were tested 
using Student's independent <-tests. General linear model (GLM) 
tests were applied to compare changes in variables between the 
two groups by adjusting for confounding factors. Pearson's and 
partial correlation coefficients were used to examine tlie relation- 
ships between variables over time. Multiple regression analyses 
were performed to identify major plasma metabolites of blood 
pressure. False Discovery Rate (FDR) corrected ^'-values were 
computed using the R package 'fdrtool'. 

Multivariate statistical analysis was performed using SIMCA-P+ 
software version 12.0 (Umetrics, Umea, Sweden). Partial 
least-squares discriminant analysis (PLS-DA) was used as the 
classification method for modeling the discrimination between 
normotensive and prehypertensive subjects by visualizing the score 
plot (S'-plot) using the first- and second-PLS components. To 
validate the model, a seven-fold validation was applied to the PLS- 
DA model, and the reliabilities of the model were rigorously 
validated by a permutation test {n = 200). The goodness of the fit 
was quantified by R^T, whereas the predictive ability was 
quantified by (fT. Generally, l^T describes how well the data 
in the training set are mathematically reproduced, and varies 
between 0 and 1 (a value of 1 indicates a model with a perfect fit). 
Models with (J^T SO. 5 are considered to have good predictive 
capabilities. 

Results 

Clinical characteristics, inflammatory and oxidative-stress 
markers, and nutrient intake of control and 
prehypertensive subjects 

The mean SBP/DBP levels in the control and prehypertensive 
groups were 107/65 and 134/84 mm Hg, respectively (Table 1). 
At baseline, the mean SBP/DBP level in the control group was 
105/65 mm Hg, compared with 134/85 mm Hg in the pre- 
hypertensive group. Subjects with prehypertension showed signif- 
icantiy higher BMI, waist-hip ratio (WHR), total- and LDL- 
cholesterol, and triglyceride than controls. The prehypertension 
group also had higher lipoprotein-associated phospholipase A2 
(LP-PLA2) activity, plasma malondialdehyde (MDA), urinary 
excretion of 8-e/)i-PGF2tx, ba-PWV, and serum interleukin 6 before 



and after adjusting for BMI, WHR, smoking, alcohol consump- 
tion, serum fipid profiles, glucose, and insulin (Table 1). The 
estimated total caloric intakes were similar in both the prehy- 
pertensive group (2,225±40 kcal/d) and the control group 
(2,194±33 kcal/d). There were no statistically significant differ- 
ences between the groups with respect to the proportion of caloric 
intake from macronutrient intake (data not shown), or the intake 
ratio of polyunsaturated/monounsaturated/saturated (P/M/S) 
fats between the control (1:0.80:0.51) and the prehypertensive 
(1:0.80:0.55) groups. There were no significant differences in the 
ratio of total energy intake to total energy expenditure between the 
groups (data not shown). 

Plasma metabolic profiling using UPLC-LTQ-Orbitrap 
mass spectrometry 

Nontargeted metabolic pattern analysis. The mass spec- 
trometry (MS) data of plasma metabolites obtained from 
normotensive and prehypertensive subjects were applied to a 
PLS-DA score plot (Figure la). The two-component PLS-DA 
score plots of the plasma metabolites showed distinct clustering 
and clear separation for each of the normotensive and prehy- 
pertensive groups. Both groups could be clearly differentiated from 
each other by the primary component t(l) or the secondary 
component t(2) based on the model with R^X{cum) and R^I{ cum 
values of 0.232 and 0.661, respectively, which indicates a good fit 
of the data. The (^3^cum) value of 0.532 provided an estimate of 
the predictive ability of the model. The PLS-DA model was 
validated using a permutation test, which indicated an R^T 
intercept value of 0.131 and a C/Tintercept value of —0.294. To 
identify the metabolites that contributed to the differentiation 
between normotensive and prehypertensive groups, tV-plots of p(l) 
and p(corr)(l) were generated using centroid scaling (Figure lb). 
The 5'-plot revealed that the metabolites with higher or lower 
p(corr) values served as the more relevant metabolites for 
discriminating between the two groups. 

Identification of plasma metabolites. The metabolites 
(variables) that play important roles in the differentiation between 
normotensive and prehypertensive groups were selected according 
to the variable important in the projection (VIP) parameter (VIP 
values >1.0 indicate a high relevance to the differences between 
the sample groups). Among 932 metabolites in plasma, 52 
metabolites were selected based on VIP values >1.0; of these, 
20 metabolites were identified and 32 were unknown. The results 
are shown in Table 2. Among the 20 identified plasma 
metabolites, the normalized peak intensities of 3 amino acids 
(leucine, phenylalanine, and tryptophan) were significantly higher 
in prehypertensive subjects than in c:ontrol subjects; however, these 
differences disappeared after adjusting for BMI, WHR, smoking, 
alcohol consumption, serum lipid profiles, glucose, and insulin. 
Also, tiiese 3 amino acids showed 0.1025, 0.1208, and 0.1000 ofq- 
value, respectively. Ten lysophosphatidylcholines (lysoPCs) con- 
taining C14:0, C16:l, C16:0, C18:2, C18:l, C18:0, C20:5, C20:4, 
C20:3, and C22:6 were significantly higher in prehypertensive 
subjects than in control subjects (9 = 5.41E-08, 5.05E-14, 6.28E- 
16, 1.27E-0.5, 4.46E-10, 1.70E-05, 6.81E-05, 6.52E-11, 2.35E-09, 
and 3.23E-09, respectively) before and after adjusting for all 
confounding variables including Lp-PLAg activity and ox-LDL. 

Blood pressure correlates with clinical and biochemical 
parameters and major plasma metabolites 

In all subjects (h = 106), SBP and DBP positively correlated with 
BMI, WHR, total- and LDL-cholesterol, triglyceride, HOMA-IR 
index, Lp-PLAj activity, ox-LDL, MDA, S-epi-PGY^^t, ba-PWV, 
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Table 1. Clinical characteristics and inflamnnatory and oxidative-stress markers. 







Normotensive Group (/7=53) 


Prehypertensive Group (/j=53) 


P 


P' 


Age (year) 


51.1 ±0.94 


51.1±1.30 


0.972 




Male/female (%) 


54.7/45.3 


54.7/45.3 


1.000 




Systolic BP (mm Hg) 


107.4±0.94 


134.1 ±0.76 


<0.001 


<0.001 


Diastolic BP (mm Hg) 


64.6 ±0.94 


85.0 ±0.56 


<0.001 


<0.001 


Body mass index (kg/m^) 


23.2±0.28 


25.3±0.37 


<0.001 




Waist/hip ratio 


0.88 ±0.01 


0.91 ±0.01 


0.001 




Cigarette smoker, n {%) 


77.4/22.6 


86.8/13.2 


0.205 




Alcohol drinker, n (%) 


30.2/69.8 


43.4/56.6 


0.159 




Total-cholesterol (mg/dL)* 


184.7 ±5.39 


208.2±4.50 


<0.001 




LDL-cholesterol (mg/dL)* 


116.1 ±4.60 


133.6±3.79 


0.001 




HDL-cholesterol (mg/dL)* 


51.7±2.05 


48.7 ±1.68 


0.277 




Triglyceride (mg/dL)* 


89.6±10.4 


133.1 ±12.4 


<0.001 




Glucose (mg/dL)* 


92.0±1.28 


95.9±1.71 


0.156 




Free fatty acid (jiEq/L)* 


468.8±32.0 


539.3±28.8 


0.050 




Insulin ((ilU/mL)* 


7.88 ±0.46 


8.92 ±0.54 


0.117 




HOMA-IR* 


1.79±0.10 


2.14±0.14 


0.067 


0.783 


hs-CRP (mg/dL)* 


1.15±0.24 


1.06 ±0.23 


0.826 


0.170 


LDL particle size (nm)* 


24.1 ±0.15 


23.5±0.15 


0.005 


0.589 


Lp-PLA2 activity (nmol/mL/min)* 


28.1 ±0.81 


33.4±0.98 


<0.001 


0.010 


Malondialdehyde (nmol/mL)* 


9.09±0.33 


11.3±0.31 


<0.001 


0.010 


Oxidized LDL (U/L)* 


38.3 ±1.30 


51.6±2.16 


<0.001 


0.051 


8-ep/-PGF2a (pg/mg creatinine)* 


1,111.4±38.8 


1,698.2±97.3 


<0.001 


< 0.001 


ba-PWV (cm/sec)* 


2,431.7±44.5 


2,99 1.3 ±58.4 


<0.001 


<0.001 


Serum IL-6 (pg/mL)" 


2.00 ±0.1 2 


3.20 ±0.1 4 


<0.001 


<0.001 


Serum IL-1 p (pg/mL)* 


0.61 ±0.10 


0.83 ±0.1 6 


0.053 


0.106 


Serum TNF-a (pg/mL)* 


6.78 ±0.73 


8.55 ±0.57 


0.046 


0.100 



Means ± S.E.*tested by logarithmic transformation. P-values derived from independent Mests. P'-values derived from independent Mests after adjusting for BIVll, WHR, 
smoking, alcohol consumption, total-cholesterol, LDL-cholesterol, HDL-cholesterol, triglyceride, glucose, free fatty acid, and insulin. H0IV1A-IR= [fasting insulin (^IU/ 
mL)xfasting glucose (mmol/L)]/22.5. hs-CRP = high sensitivity C-reactive protein. ba-PWV = brachlal-pulse wave velocity. 
doi:1 0.1 371 /journa!.pone.0096735.t001 



IL-6, lysoPCs (14:0, 16:1, 16:0, 18:2, 18:1, 18:0, 20:5, 20:4, 20:3, 
and 22:6) (aU Pvalues <0.01), and IL-ip (SEP, P=0.040; DBF, 
P= 0.020), and negatively correlated with LDL particle size (all P 
values <0.01). Systolic blood pressure also positively correlated 
with tryptophan {P= 0.01.5). Based on these results, we performed 
a multiple-regression analysis to determine independent predictors 
of SBP and DBF. Age, sex, BMI, WHR, total cholesterol, HOMA- 
IR index, Lp-PLA2 activity, leucine, phenylalanine, tryptophan, 
and lysoPCs (14:0, 16:1, 16:0, 18:2, 18:1, 18:0, 20:5, 20:4, 20:3, 
and 22:6) were tested. LysoPC (16:0) emerged as an independent 
predictor of SBP (standardized ;8 = 0.399, P= 0.046), as did BMI 
(standardized /? = 0.225, P= 0.025). LysoPC (16:0) also was an 
independent predictor of DBP (standardized P = 0.453, P= 0.025). 

LysoPC (16:0) correlates with ox-LDL, Lp-PLAz activity, 8- 
ep;-PGF2a, IL-6, IL-ip, TNF-ot, and other lysoPCs 

In prehypertensive subjects, the levels of lysoPC (16:0) positively 
and significantly correlated with ox-LDL, Lp-PLA2 activity, 8-epi- 
PGFga, and IL-6 before and after adjusting for BMI, WHR, 
smoking, alcohol consumption, and LDL cholesterol (Figure 2). 
However, these associations were not found in normotensive 
subjects. In both prehypertensive and control subjects, the levels of 
lysoPC (16:0) positively correlated with other lysoPCs (14:0, 16:1, 



18:2, 18:1, 18:0, 20:5, 20:4, 20:3, and 22:6) and ba-PWV (aU P 
values <0.014). Serum IL-6, IL-1|3, and TNF-a positively 
correlated with urinary H-epi-PGF^a (r= 0.262, P= 0.007; 
)■= 0.281, P= 0.042; r= 0.373, P= 0.006, respectively) and plasma 
MDA ()■= 0.449, P<0.001; ) = 0.279, P= 0.043; r= 0.277, 
P = 0.044, respectively). 

Discussion 

This study identified prehypertension-associated alterations in 
lysoPC levels and amino acid metabolism. Significant differences 
in the metabolic profiles were found between prehypertensive and 
normotensive individuals, including lysoPCs (14:0, 16:1, 16:0, 
18:2, 18:1, 18:0, 20:5, 20:4, 20:3, and 22:6) and amino acids 
(leucine, phenylalanine, and tryptophan). Among the 10 lysoPCs, 
the majority contained long-chain acyl groups (i.e., C&16) [12]. 
This result is in agreement with previous studies on the elfects of 
long-chain lysoPC (CS16) in vasodilation impairment [13], and 
the identification of higher levels of lysoPCs with long-chain acyl 
groups in the plasma of spontaneously hypertensive rats.** This 
study identified lysoPC (16:0) (VIP value of 17.173) as the most 
important plasma metabolite for evaluating the difference between 
prehypertensive and normotensive individuals. LysoPC (16:0) was 
positively and independently associated with DBP and SBP. This 
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Figure 1. Partial least-squares discriminant analysis models. (A) Score plots classifying subjects as normotensive (filled square) or 

prehypertensive (filled triangle). (B) Score plots for covariance [p] and reliability correlation [p(corr)]. 
doi:1 0.1 371 /journal.pone.0096735.g001 



observation could support previous work on the role of increased 
lysoPC (16:0) in arterial stiffness [14]. 

LysoPC constitutes only 1-5% of the total PC content of non- 
ox-LDL; however, as much as 40—50% of the PC contained within 
LDL is converted to lysoPC during LDL oxidation [15]. A 



saturated or monounsaturated fatty acid predominates in the sn- 1 
position of the phospholipid [16]. The generation of free radicals 
as a result of oxidative stress can activate PLA2, which catalyzes 
the hydrolysis of the ester bond at the sn-2 position of 
phospholipids. Stafforini et al. [16] showed that the secreted form 
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of Lp-PLAg released F2-isoprostanes (the end-products of lipid 
oxidation) from the sn-2 position of PC with high affmity. Kono et 
al. [1 7] reported that intracellular type II Lp-PLA2, which shares 
homology with the plasma enzyme Lp-PLA2, was involved in the 
metabolism of esterified 8-ifo-PGF2o,. We also measured prehy- 
pertension-associated elevations in plasma Lp-PLA2 activity, 
plasma ox-LDL, lysoPCs, and urinary 8-e/ii-PGF2oc These results 
are in agreement with a previous report of higher mean levels of 
Lp-PLA2 ^i^d 8-epi-?GF2a in subjects with high ox-LDL than those 
in subjects with low ox-LDL [18], and a positive association 
between plasma ox-LDL and Lp-PLA2 activity in metabolic 
syndrome [19], 

In prehypertensive individuals, a strongly positive correlation of 
lysoPC (16:0) with ox-LDL, Lp-PLA2 activity, and 8-epi-PGF2a (a 
sensitive marker for oxidative stress) [20,21] could indicate 
increased production of oxidative stress from Lp-PLA2-catalyzed 
PC hydrolysis during increased LDL oxidation. This result is 
consistent with a previous report of a correlation between 
increased Lp-PLA2 activity and increased levels of lysoPC, ox- 
LDL, and cytokine in hypercholesterolemic minipigs [22]. A close 
correlation between the circulating lysoPC (16:0) levels and 
urinary 8-ej()f-PGF2oi also was observed in middle-aged nonobese 
men [14]. However, this association was not found in normoten- 
sive subjects, possibly due to low levels of LDL cholesterol, ox- 
LDL, and Lp-PLA2 activity. 



LysoPC, a major atherogenic component of oxLDL, was 
elevated in the plasma of atherosclerotic patients [23], and was 
an important biomarker of coronary artery disease (CAD) [24]. 
Lp-PLA2 levels were observed to correlate with three lysoPCs 
(16:0, 18:0, and 18:1) in carotid plaques before and after adjusting 
for age, sex, hypertension, CAD, and statin usage [10]. A close 
association between circulating lysoPC (16:0) and the ba-PWV 
index of central arterial stiffness [25,26], and serum IL-6, was 
reported in human subjects [14]. Production of lysoPC (16:0) can 
stimulate the expression of adhesion molecules and the release of 
cytokines in endothelial cells [26]. IL-6 was induced by lysoPC 
(16:0) treatment in human umbilical vein endothelial cells. LysoPC 
(16:0) induced higher release of arachidonic acid mediated via 
cytosolic PLA2 from human coronary artery smooth muscle cells, 
compared to that for treatment with lysoPC (14:0) or unsaturated 
lysoPC [27]. LysoPC has been also shown enhance the release of 
IL-6 from peripheral blood mononuclear cells [28], regulated on 
activation normal T-cell expressed and secreted from microvas- 
cular endothelial cells [29], and TNF-a [30] as well as IL-lfi from 
human monocytes [31]. In addition, previous epidemiologic 
studies have linked higher plasma concentrations of inflammatory 
markers - IL-6 [32,33], IL-1(3 [34], and TNF-a [35] - to increased 
SBP and DBP or hypertensive status. Consistent with these 
reports, the current study showed significant increases in ba-PWV, 
serum IL-6 levels, and lipid peroxides including 8-epi-J'G¥.2oi and 
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MDA in the prehypi^rtcnsivc group, and positive correlations 
among lysoPC (16:0), IL-6, 8-epi-FGF2^, and MDA. 

LysoPCs represent 5-20% of the total plasma phospholipids, 
and are formed by the action of lecithin cholesterol acyltransferase 
(LCAT) in plasma [36]. Human LCAT releases lysoPC 20:4 and 
22:6 from the sn-l position of PC [37]. Up to 80% of the lysoPC in 
plasma is found in the non-lipoprotein fraction, in which albumin 
is considered as the main hpid-binding protein [36,38]. Unsatu- 
rated lysoPGs are associated primarily with albumin rather than 
lipoproteins. We observed a positive relationship between the 
levels of lysoPC (16:0) and those of other lysoPCs (14:0, 16:1, 18:2, 
18:1, 18:0, 20:5, 20:4, 20:3, and 22:6), which could reflect an 
alternative source of lysoPC (16:0) production in addition to ox- 
LDL. LysoPC is formed by PLAg-induced hydrolysis or oxidation 
of PC in LDL and cell membranes [23]. LysoPCs (14:0, 16:0, 18:0, 
and 18:1) and branched-chain or aromatic amino acids were 
elevated in overweight/obese, insulin-resistant subjects compared 
with the levels in lean control subjects [39-42]. In the current 
study, differences in the levels of leucine, phenylalanine, and 
tryptophan between control and prehypertensive subjects disap- 
peared after adjusting for BMI, WHR, smoking, alcohol 
consumption, serum lipid profiles, glucose, and insulin. However, 
differences in the levels of lysoPCs in the two groups remain highly 
significant. This result suggests an important role of lysoPCs in 
prehypertension, independendy from confounding variables. 

Our results have similar limitations as all cross-sectional and 
observational studies. We evaluated associations rather than 
prospective prediction; thus, the causal relationships among the 
identified metabolites and the exact mechanisms of the prehy- 
pertensive changes are unknown. A large number of markers were 
detected using UPLC-LTQ;Orbitrap MS, but most remain 
unidentified. Unlike the large databases for GC-MS, the 
endogenous biomolecule databases for LC-MS-based metabolo- 
mics research have not yet been constructed [43] . Our study used 
UPLC-LTQ;Orbitrap MS-based metabolomics and multivariate 
data analyses to identify a cluster of prehypertension-associated 
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